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Abstract 

The longicorn beetles belonging to the genus MesechthistatusBreuning, 1950 (Coleoptera: Cerambycidae) cannot fly since 
their hindwings are atrophied. This slows down gene flow between local populations. Currendy, it is considered that the 
genus contains four endemic species from the eastern Honshu Is., Japan, M. binodosus, M.furciferus, M. taniguchii and M. 
fujisanus, that are distributed parapatrically. Sequence analyses of the cytochrome oxidase subunit I gene suggests that 
lineages of mitochondrial haplotypes split approximately in the same era. However, this result is not consistent with the 
monophyly of morphological species. The estimated evolutionary rate of the COI gene in other insects suggests that 
mitochondrial haplotypes of Mesechthistatus differentiated at the end of the Pliocene epoch during the Tertiary era. 
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Introduction 

Situated at the eastern tip of the Eurasian plate, the 
Japanese Archipelago forms an island arc with a total 
length of 3000 km (Figure 1). The landmass was separ- 
ated from the Eurasian continent between 20 and 16 mil- 
lion years ago (mya) during the Miocene epoch of the ter- 
tiary period of the Cenozoic era and became the basic 
elements of the present day archipelago (Maruyama et al. 
1997; Otofuji et al. 1985, 1996). Subsequently this land 
sank, and the land started to rise in the west, resulting in 
the formation of the archipelago around 5 million years 
ago in the Pliocene epoch (Maruyama et al. 1997). More 
recent studies (Kitamura et al. 2001; Kitamura and 
Kimoto 2006) reported that the Japanese Archipelago 
was connected to the East Asian continent between 3.5 
and 1.7 mya except for a few short periods of isolation 
lasting 10,000 to 20,000 years. After 1.7 mya, in the 
Quaternary era, it is considered that unstable land 
bridges were constructed several times between the 
Japanese Archipelago and the Continent at glacial max- 
ima (Kawamura 2007; Shinohara et al. 2004; Suzuki et 
al. 1997, 2003; Tsuchiya et al. 2000). These caused geo- 
graphic links and separation from the continental biome, 
and subsequent colonization of the archipelago com- 
pleted the speciation process between the early and late 
colonizers in the archipelago. In addition, the rising land 
level in the period after the Quaternary era in the central 
region of Honshu resulted in complicated topography. 
For animals of low mobility, high mountains and other 
geographical features became barriers to migration, and 
it is thought that separation and speciation occurred at 
regional population levels. 

These events in the geological history formed a diversity 
of endemic species in the Japanese Archipelago. For ex- 
ample, an extensive taxonomic survey described 755 spe- 
cies of Cerambycidae, of which 463 are endemic to 
Japan (Ohbayashi and Niisato 2007). Investigations to 
the subfamily level gave the following figures (the number 
of endemic species/the number of total species in the 
Japanese Archipelago): Parandrinae; 1/2, Prioninae; 
4/13, Spondylidinae; 4/15, Lepturinae; 118/157, Necy- 
dalinae; 7/10, Cerambycinae; 97/198, Lamiinae; 232/ 
360 with notably high numbers of endemic species in 
Lepturinae and Lamiinae. 

The genus MesechthistatusBreuning, 1 950 discussed in this 
paper belongs to the tribe Phrissomini Thomson, 1860 
and the subfamily Lamiinae Latreille, 1825, and contains 
four species based on morphology and distribution pat- 
terns: M. binodosus (W aterhouse), M. furciferus (Bates), M. 
taniguchii (Seki, 1944), and M. fujisanus Hayashi, 1957 
(Figure 2). Mesechthistatus jamahoi (Mitono, 1943), previ- 
ously described in Taiwan (Mitono 1 943), also belongs to 
this genus. However, with no subsequent records since 
the time of original entry, the existence of this species is 
currendy in doubt (Hasegawa 2007). Mesechthistatus is a 
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genus endemic to Japan, and is found only from Honshu 
and Sado Island. The above-mentioned four species have 
parapatric distribution (Figure 1). The atrophied hind- 
wings common to all species render them incapable of 
flight, and thus their mobility must be limited. As a res- 
ult, morphological disparities have been found, with two 
subspecies for each of the species described as M. 
binodosus and M. furciferus (Hayashi 1951, 1955). 

In this paper, to investigate interspecific differences of 
Mesechthistatus and to estimate the time of colonization in 
Japanese Archipelago, molecular phylogenetic analysis 
was carried out using partial sequences from the mito- 
chondrial cytochrome oxidase subunit I gene from four 
different species of Mesechthistatus. 

Materials and Methods 

Taxon sampling 

The specimens of Mesechthistatus spp. analyzed in this 
study are listed in Table 1, and the localities where they 
were collected are shown in Figure 3. The beedes were 
immediately fixed in 95-99.5% ethanol and preserved in 
the same solution until use. A single individual from each 
locality was used for DNA extraction. 

Extraction, PCR amplification and sequence 
analysis of DNA 

Total DNA was extracted from a mixture of cephalic and 
thoracic muscles by using a GenElute™ Mammalian Ge- 
nomic DNA Miniprep Kit (Sigma-Aldrich Inc., 
www.sigmaaldrich.com). Each DNA sample was dis- 
solved in 200 |x"l elution buffer. A fragment of DNA en- 
coding mitochondrial cytochrome oxidase subunit I was 
amplified from the total DNA solution using PCR with a 
primer pair as follows: KobCI1.2 (5 -TAA GAA GAA 
TTG TAG AAA ATG G-3') and YhzCI2.2 (5'-TGT 
AGC GAT TTC TAA AAA AAGG-3'). PCR was car- 
ried out in a 25 |x"l reaction mixture containing 1 X 
PCR buffer for KOD -Plus- (Toyobo Biologies Inc., 
www.toyobobiologics.com), 0.2 mM of each dNTP, 2 
mM of MgS04, 0.5 unit of KOD -Plus- DNA Poly- 
merase (Toyobo), 0.3 |x"M of each primer and 1 u,"l of 
template DNA solution. The amplification protocol was 
25 cycles of denaturation at 94°C for 15 sec, annealing at 
50°C for 30 sec, and extension at 68°C for 40 sec in a 
GeneAmp® PCR System 9700 (Applied Biosystems, 
www.appliedbiosystems.com). The PCR product was 
purified using by a GenElute™ PCR Clean-Up Kit 
(Sigma-Aldrich Inc.). Direct sequencing of the COI DNA 
fragment was performed by using a BigDye® Terminat- 
or Cycle Sequencing Kit (Applied Biosystems) with 
primers KobCI1.2 and Yhz2.2. A partial sequence of the 
COI gene was determined by a ABI PRISM® 310 Gen- 
etic Analyzer or ABI PRISM® 3100 Genetic Analyzer 
(Applied Biosystems). 
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Figure I.Japan showing the distribution areas of four Mesechthistatus species. 



Phylogenetic analysis 

The sequence alignments were performed by using 
CLUSTAL W version 1.83 (Thompson et al. 1994), and 
no indels (insertions/deletions) were required for the 
alignments of COI sequences. A transition/transversion 
ratio R was obtained with PAUP* version 4. Ob 10 
(Swofford 2002). 

Phylogenetic tree topology was constructed according to 
the maximum likelihood (ML) method in PAUP* version 
4.0b 10. The ML model was selected by hierarchical like- 
lihood ratio tests (hLRTs) and akaike information cri- 
terion (AIC) using MrModeltest version 2.2 (Nylander 



2004) with PAUP* version 4.0b 10, and heuristic searches 
were performed using the GTR+I+G model of the sub- 
stitutions with TBR branch-swapping algorithm. The 
starting tree was obtained via stepwise addition, and the 
starting branch lengths obtained using Rogers-Swofford 
approximation method. The bootstrap test was executed 
on 200 replicates using GARLI version 0.951 (Zwickl 
2006). 

Estimation of divergence time 

To estimate divergence time, PATHd8 was used (Britton 
et al. 2007). The PATHd8 is a program for phylogenetic 
dating without a molecular clock. First, a ML tree 
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Figure 2. Photographs of Mesechthistatus spp. I , M. binodosus; 2, M. furciferus; 3, M. taniguchii; 4, 
M. fujisanus 



the upper fix age and 0.2 mya was used as the lower fix 
age at the node of "BIN 8" and "BIN 12". 

Results and Discussion 

Genetic data 

A partial sequence of 1144 base pairs of mitochondrial 
COI gene was obtained from Mesechthistatus spp. In this 
study, throughout the COI sequences used, neither dele- 
tions nor insertions were found in multiple alignment. Al- 
together 236 sites were variable including 33 at the first 
position, 4 at the second position and 189 at the third po- 
sition. The transition/transversion ratio R was 5.103. 
The base frequencies were nearly constant among four 



topology was prepared using PAUP* version 4.0b 10 with 
the GTR+I+G model. Next, branch lengths were ob- 
tained using an estbranches program that included the 
multidistribute program package (Thorn et al. 1 998) with 
PAML (Yang 1997). Thus, we prepared input file (tree 
with branch lengths in Newick format) for PATHd8. 

One constraint was used on one node. Mesechthistatus 
binodosus is distributed in Honshu and Sado Is. (Figure 1). 
Oshima (1990) estimated that a channel separated Sado 
Island and Honshu between 0.8 and 0.2 Mya. The ob- 
tained likelihood branch length between "BIN 12" col- 
lected from Sado Island and "BIN 8" collected from 
Honshu, are the shortest. Therefore, 0.8 mya was used as 
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Table I. The specimens analyzed in this study. 



Scientific Name (by 
morphology) 


Isolate 
Code 


No. in 
Map 


Locality 


DDBJ/EMBL/GenBank 
Accession number 


Mesechthistatus binodosus 




BIN 


1 


Aiuchi, Shiura-mura, Aomori 


AB27822I 




BIN 


2 


Anmonnotaki, Nishimeya-mura, 
Aomori 


AB278235 




BIN 


3 


Genbeidaira, Niisato-mura, Iwate 


AB278240 




BIN 


4 


Mt. Moriyoshi, Ani-machi, Akita 


AB278247 




BIN 


5 


Mt. Moriyoshi, Ani-machi, Akita 


AB278248 




BIN 


6 


Shinzan, Oga-shi, Akita 


AB278249 




BIN 


7 


Mt. Hinoto, Chokai-machi, Akita 


AB278252 




BIN 


8 


Nabekoshi Pass, Obanazawa-shi, 
Yamagata 


AB278253 




BIN 


9 


Futakuchi Pass, Yamagata-shi, 
Yamagata 


AB278258 




BIN 


10 


Kotamagawa, Oguni-machi, 
Yamagata 


AB278259 




BIN 


1 1 


Usokurasawa, Minakami-machi, 
Gunma 


AB278263 




BIN 


12 


Mt. Donden, Aikawa-machi, Sado 
Is., Niigata 


AB278265 




BIN 


13 


Mt. Taho, Iwamuro-mura, Niigata 


AB278268 




BIN 


14 


Sasagamine, Myokokogen-machi, 
Niigata 


AB278276 




BIN 


15 


Sasagamine, Myokokogen-machi, 
Niigata 


AB278277 




BIN 


16 


Mt. Narou, Suzaka-shi, Nagano 


AB278297 




BIN 


17 


Takeshi Pass, Takeshi-mura, 
Nagano 


AB278308 




BIN 


18 


Sano, Hakuba-mura, Nagano 


AB2783I2 




BIN 


19 


Taira, Omachi-shi, Nagano 


AB2783I4 


M. furciferus 






FUR 


20 


Mt. Oitsurushi, Nyuzen-machi, 
Toyama 


AB278325 




FUR 


21 


Mt. lou, Fukumitsu-machi, 
Toyama 


AB27833 1 




FUR 


22 


Mt. Horyu, Wajima-shi, Ishikawa 


AB278333 




FUR 


23 


Mt. Kariyasu, Awara-shi, Fukui 


AB278337 




FUR 


24 


Mt. Hoonji, Katsuyama-shi, Fukui 


AB278338 




FUR 


25 


Mt. Jo, Echizen-cho, Fukui 


AB278346 




FUR 


26 


Mushio, Otari-mura, Nagano 


AB278359 




FUR 


27 


Oyasawa right bank, Otari-mura, 
Nagano 


AB278363 




FUR 


28 


Oyasawa right bank, Otari-mura, 
Nagano 


AB278364 




FUR 


29 


Odairahara, Hotaka-machi, 
Nagano 


AB278373 




FUR 


30 


Hiwadakogen, Takane-mura, Gifu 


AB278374 




FUR 


31 


Oshirakawa, Shirakawa-mura, 
Gifu 


AB278377 



Journal of Insect Science | www.insectscience.org 



5 



Journal of Insect Science: Vol. 8 | Article 70 



Nakamine et al. 



Table I (con't.) 



Scientific Name (by 
morphology) 


Isolate 
Code 


No. in 
Map 


Locality 


DDBj/EMBL/GenBank 
Accession number 




FUR 


32 


Kanmuriyama Pass, 
Fujihashi-mura, Gifu 


AB27838I 




FUR 


33 


Mt. Oiks, Fujiwara-cho, MiG 


AB278384 




FUR 


34 


Kuroko Pass, Makino-cho, Shiga 


AB278388 




FUR 


35 


Asyu, Miyama-cho, Kyoto 


AB278398 




FUR 


36 


Mt. Choro, Wachi-cho, Kyoto 


AB278399 




FUR 


37 


Mt. Oe, Oe-cho, Kyoto 


AB278400 




FUR 


38 


Uchiyama, Oomiya-cho, Kyoto 


AB27840I 




FUR 


39 


Nasa, Toyooka-shi, Hyogo 


AB278403 




FUR 


40 


Itoikeikoku, Wadayama-cho, 
Hyogo 


AB278406 


M. taniguchii 






TAN 


41 


Maruyama rindo, Masuho-cho, 
Yamanashi 


AB2784I 1 




TAN 


42 


Yadehara rindo, Kawakami-mura, 
Nagano 


AB2784I3 




TAN 


43 


Mt. Torikura, Oshika-mura, 
Nagano 


AB2784I6 




TAN 


44 


Sengendaru, Takane-mura, Gifu 


AB278427 




TAN 


45 


Kuroisawa rindo, 
Nakatsugawa-shi, Gifu 


AB278428 




TAN 


46 


Mt. Yanbushi, Shizuoka-shi, 
Shizuoka 


AB278430 


M. fujisanus 






FUJ 


47 


Mt. Haruna, Agatsuma-machi, 
Gunma 


AB278434 




FUJ 


48 


Dodaira, Kiyokawa-mura, 
Kanagawa 


AB278437 




FUJ 


49 


Otomekogen, Koufu-shi, 
Yamanashi 


AB27844I 




FUJ 


50 


Takinoirisawa, Sanada-machi, 
Nagano 


AB278445 




FUJ 


51 


Nishiusuzuka, Fujinomiya-shi, 
Shizuoka 


AB278450 


Out group 
Parechthistatus gibber 






GIB 


52 


Mt. Oto, Kozagawa-cho, 
Wakayama 


AB278487 


Hayashiechthistatus inexpectus 




INE 


53 


Yodogawa, Yaku Is. 
Kamiyaku-cho, Kagoshima 


AB27855 1 



Mesechthistatus spp. as follows: A, 0.309971; C, 0.163393; 
G, 0.149304 and T, 0.377331. These data suggest an A- 
T bias. 

Molecular phylogenetic analysis 

Figure 4 shows the maximum likelihood phylogenetic 
tree. Basically four lineages of ancestry were recognized, 



as marked A, B, C and D. A schematic diagram for re- 
latedness of the 4 lineages is shown in Figure 5. 

Lineage A: Lineage A contained a single species, M. 
binodosus, showing the monophyletic origin of the M. 
binodosus COI haplotype (Figure 4). However, morpholo- 
gical characteristics and the fact that the location of col- 
lection was within that of the Al. furciferus distribution 
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Figure 3. The localities where specimens of four Mesechtfi/status species were collected. Open circle, locality of M. binodosus speci- 
mens; closed circle, locality of M. furciferus specimens; open square, locality of M. fujisanus specimens; closed square, locality of M. 
taniguchii specimens; closed star, localities of Parechthistatus and Hayashiechthistatus used as outgroup specimens for phylogenetic ana- 
lysis. Locality numbers correspond to the number in Table I and the phylogenetic trees in Figure 4. 



area imply that M. furciferus (FUR 28) should be included 
in linage A. As a specimen of Al. furciferus was collected in 
the vicinity of the M. binodosus distributional boundary, it 
is likely that introgressive hybridization (e. g. Nagata et 
al. 2007; Shimizu and Ueshima 2000; Sota et al. 2001) 
took place from Al. binodosus to Al. furciferus sometime in 
the past, thus giving "FUR 28" the Al. binodosus haplo- 
type. Given that the "FUR 27" that was collected at the 
same location as "FUR 28" has the same COI haplotype 
as Al. furciferus, it is likely that individuals that underwent 
introgressive hybridization now exist together as a single 
group in one location. 



Lineage B: This lineage was not monophyletic but com- 
plex, including 3 species: Al. furciferus, Al. tanguchii, and Al. 
fujisanus. 

B-l contains a single species, Al. furciferus and shows a 
unique distribution pattern (Figure 5). Lineage D is dis- 
tributed in the gap zone between the B-l populations 
since it has another ancestral mitochondrial haplotype. 
Lineage D is considered to have diverged at an earlier 
time than the lineage B, and it is possible that the mito- 
chondrial haplotype of lineage D expanded its distribu- 
tion into the area of lineage B distribution. Lineage B di- 
verged into several sublineages thereafter. 
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Figure 4. Maximum likelihood tree based on the mitochondrial COI gene sequences of the Mesechthistatus 
spp. The bootstrap value was indicated at each node (when > 50%). The numbered specimen identity cor- 
responds to the number for locality in Table I, Figure 3 and Figure 5. Asterisk indicates introgressive 
specimens. 



M. furciferus, M. taniguchii, and M. fujisanus form the B-2 
sublineage. The B-2-a sublineage is made of populations 
found in the vicinity of the western limit of the M. 



furciferus distribution. The B-2-b sublineage is essentially 
made of M. taniguchii showing monophyletic origin, 
though an individual "FUJ 48" was recognized 
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Figure 5. A simplified tree showing relationship between major lineage and geographic distribution in the study. The letter beside 
each lineage line indicates correspondence to the letter in Figure 4. The bootstrap value was indicated at each node (when > 50%). 
Open circle, locality of M. binodosus specimens; closed circle, locality of M. furciferus specimens; open square, locality of M. fujisanus 
specimens; closed square, locality of M. taniguchii specimens. 



morphologically as Al. fujisanus. Until now, neither Al. 
taniguchii nor Al. fujisanus has been found to possess a mi- 
tochondrial haplotype closely related to the haplotype 
found in this sample. This suggests introgressive hybridiz- 
ation between Al. taniguchii and Al. fujisanus at some time 
in the past. Sublineage B-2-c includes M. fujisanus, and 
the genetic disparity among populations was small. Poss- 
ible explanations for this include the emergence of a bot- 
tleneck and rapid distribution expansion at some time in 
the past. 

Within the B-2 sublineage, there is a contiguity in the dis- 
tribution areas of B-2-b (Al. taniguchii), and B-2-c (Al. fujis- 
anus). However, the distribution area of the B-2-a popula- 
tion occurs at the western border of Al. furciferus, and it is 
isolated from the other lineage by a considerable dis- 
tance. Mesechthistatus has not been found in the southwest 
region of Honshu, Shikoku and Kyushu, although the 



genus Parechthutatm, a member of the tribe that Mesech- 
thistatus belongs to, inhabit that region. It is possible that 
the groups with the mitochondrial haplotype found in the 
B-2-a became fragmented due to the distribution expan- 
sion of Parechthistatus. 

Lineage C: Al. furciferus and Al. binodosus form the lineage 
C from an extremely small area in the northern part of 
the Chubu region (Figure 5). It is considered that this lin- 
eage originally possessed the M. furciferus mitochondrial 
haplotype, though "BIN 15" and "BIN 19" are included 
in this lineage. As this population is found in the vicinity 
of the distribution boundary, it is likely that Al. furciferus 
mitochondrial DNA was introgressed into Al. binodosus. 

Lineage D: Al. furciferus forms a single lineage D, distrib- 
uted around Lake Biwa (Figure 5). 
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The divergence date of Mesechthistatus spp 

For molecular phylogenetic analysis, estimation of the 
evolutionary rate of sequence divergence is of utmost im- 
portance as a molecular clock. The estimation of the final 
divergence time based on this rate is a major goal. 
However, there is little knowledge of geological events 
and a lack of fossil records for calibration of the chrono- 
logical relation between phyletic lineages and genetic dis- 
tances in the Mesechthistatus COI gene. Therefore, it is dif- 
ficult to exacdy estimate the rate of evolution of the COI 
gene in this study. However, the calibration of the evolu- 
tionary rate has been estimated based on COI gene in 
numerous other insects, giving an estimate ranging from 
1.5% per 1 million years in Tetraopes: Cerambycidae 
(Farrell 2001) and Crematogaster. Formicidae (Quek et al. 
2004), 1.6% in Plateumaris: Chrysomelidae (Sota and 
Hayashi 2007) to 2.3% in Heliconius: Nymphalidae 
(Brower 1994). The date of radiation in Mesechthistatus 
was estimated by using these values. 

Figure 4 depicts four major lineages, A to D. The aver- 
age value for genetic divergence between the four lin- 
eages was 4.53 + 0.45 (mean + SD) %. Application of 
the evolutionary rate of the COI gene (1.5% to 2.3%) 
dates the divergence timing of the major four lineages 
between 3.02 and 1.96 mya. 

The divergence times were also estimated using the 
PATHd8 program that is not based on a molecular 
clock. The divergence time at the root node was between 
5.44 and 1.36 mya, in the constraints that Sado Island 
was formed between 0.8 mya, as the upper fix age, and 
0.2 mya, as the lower fix age, which corresponds to the 
distance for node of "BIN 8" and "BIN 12". This result 
does not contradict the estimated times indicated above. 

It is probably appropriate to assume that the ancestor of 
Mesechthistatus colonized from the East Asian continent, 
because Parechthistatus chinensisY^rtunmg, 1 942 that was re- 
lated to Mesechthistatus, has been recorded in the Shaanxi 
Province of the People's Republic of China (Breuning 
1942). Therefore, it is considered that colonization by 
Mesechthistatus of the Japanese Archipelago occurred in 
the late Pliocene epoch. The Japanese Archipelago was 
connected to the East Asian continent between 3.5 and 
1.7 mya (Kitamura et al. 2001; Kitamura and Kimoto 
2006). 

Sota and Hayashi (2007) reported colonization history of 
Plateumaris Thomson, 1859 leaf beedes in Japan. They es- 
timated the timing of colonization by using the Bayesian 
approach based on fossil records calibrations. They con- 
cluded that Plateumaris constricticollis Jacoby), a Japanese 
endemic species, was formed before the late Pliocene. 

Therefore, these results and other reports suggest that 
Mesechthistatus underwent radiation at the end of the Plio- 
cene epoch. In order to further investigate the evolution 



of speciation in Mesechthistatus, it will be necessary to con- 
duct a comprehensive molecular phylogenetic analysis 
that includes nuclear DNA. 

Editor's note 

Paper copies of this article will be deposited in the follow- 
ing libraries. Senckenberg Library, Frankfurt Germany; 
National Museum of Natural History, Paris, France; 
Field Museum of Natural History, Chicago, Illinois USA; 
the University of Wisconsin, Madison, USA; the 
University of Arizona, Tucson, Arizona USA; Smithsoni- 
an Institution Libraries, Washington D.C. USA; The 
Linnean Society, London, England. 
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